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ABSTRACT 

Space weathering is the generic term used for processes that modify the optical 
properties of surfaces of atmosphereless rocky bodies under exposure to the space 
environment. The general agreement about the relevance of the effects of space 
weathering on the spectral properties of S-complex asteroids fails when some basic 
quantitative estimates are attempted. In particular, there is severe disagreement 
regarding the typical timescales for significant spectral reddening to occur, ranging 
from 1 Myr to 1 Gyr. 

Generally speaking, the spectral reddening of an individual object can be considered 
as the sum of three terms, one (which is relevant for statistical analyses) depending 
on the exposure of the object to space weathering during its lifetime, a second one 
due to the original surface composition, and a third one (a noise term) due to the 
combination of poorly constrained effects (e.g., structure and texture of the surface). 
The surface of an asteroid is usually covered by regolith, and its presence and 
properties presumably play a critical role in the weathering processes. In this paper 
we discuss the role played by collisional evolution in affecting the spectral properties 
of asteroids and refreshing the surfaces due to the formation of ejecta, and the 
necessity of a simultaneous modeling of collisions and weathering processes. We 
introduce a new idea, based on the possibility of a sort of saturation of the refreshing 
process whenever a massive reaccumulation of the impact ejecta takes place. In this 
case, a dependence of the overall reddening on the asteroid size should naturally come 
out. We show that this conclusion is indeed supported by available main belt asteroid 
spectroscopic data. 



Key words: minor planets, asteroids: general; meteorites, meteors, meteoroids 



1 INTRODUCTION 

The spectral reddening of asteroidal surfaces is of 
fundamental importance to understanding the spectral 
properties derived from remote observations of asteroids 
and for making a reliable comparison with the correspond- 
ing laboratory meteorite data. The process of spectral 
reddening for S-complex asteroids has been r ecognized for 
a lon g time (|Chapman and Salisbury! 1 19731 : iGaffev et al.l 
119931 ). and it has been attributed to the alteration of 
surface optical properties under exposure to the space 
environment. More rec ently, several contributions from 
labor atory experiments l|Hiroi et al.l 1200 ll : IStrazzulla et al.l 
I2OO5I ). in addition to dedicated asteroid observations and 



modeling, have disclos ed several intr i guing a spects of the 
reddening processes Jjedicke et al.l |2004| : iMarchi et al.l 
l2006al : iPaolicchi et al. I l2007l ). It has also been shown that, 
to some extent, the overall results obtained from S-complex 
asteroids may play a non- negligible role also for other main 
taxonomic complexes a nd spectral types (jLazzarin et al.l 
I2OO6I : iMarchi et al]|201Gl ) . 

On the basis of a statistical analysis of a large sample 
of visi ble spectra of S-complex asteroids IMarchi et al.l 
((iooii) (hereinafter PI), suggested that the solar wind is 
the dominant cause of the reddening. This conclusion was 
based on the significant dependence of the spectral slope on 
the exposure (E) to the solar wind, namely the integrated 
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ion flux that an asteroid received from the Sun during its 
past evolution. This analysis is based on a large sample 
of both near-Earth asteroids (NEAs) and main belt aster- 
oids (MB As). In a series of papers (|Marchi et al.l l2006bl : 
iPaolicchi et"al]|2007l ). additional details of the process were 
determined, and an estimate of the reddening timescale 
was given. It was found that 80% of the slope reddening is 
reached after about 2 00 Myr at 1 AU; or 800 Myr at 2 AU 
(jPaolicchi et al.ll2007l ). In general, the derived correlation 
between the slope and the asteroid age, the latter estimated 
by the coUisional lifetime, requires that the time derivative 
exhibit an exponential decay: a fast initial reddening 
(leading to almost 50% of the final reddening in a few 
tens Myr) is followed by a slower increase, asymptotically 
pointing to saturation. Note, however, that the estimated 
reddening timescale due to heavy ion bo mbardment is of 
the order of 0.01 - 1 Myr at 1 AU jStrazzu lla ct al. 2005). 
A recent analysis b ased on S-complex MBA families 
jVernazza et al.l |2009| , hereinafter P2) claims that the 
solar wind is the dominant cause of the reddening. This 
conclusion was derived on the basis of the fast reddening 
timescale (~ 1 Myr) observed for two young asteroid 
families. Moreover, a slower residual reddening is present 
in their data for older families. Despite P2 agrees on the 
what already found by PI about the solar wind-driven 
space weathering, the derived timescales differs by orders 
of magnitude. 

To make the s t ory m ore complicated, a recent analysis 
l|WiUman et all I2OI0I ') suggests that space weathering 
timescales may be very large (of the order of 10^ — 10"^ Myr 
for MB As) , claiming agreement with a different series of lab- 
oratory experiments that assume micromete orite bombard- 
ment is the dominant weathering process ijWillman et al.l 
I2OI0I . and references therein). While the timescales found 
by Willman et al. are partially in agreement with some PI 
results, the physical process that is the claimed cause of the 

reddening is completely different. 

A recent comprehensive review by iGaffevI ^oM) 

points 

out that the consequences of space weathering may vary 
strongly among bodies, due, possibly, both to different 
physical processes and to different properties of the affected 
surfaces. For instance, it has been found that the space 
weathering of lunar regolith is likely due to the formation 
of nanophase iron particles (|Pieters et al.l l2000 l) produced 
by the vaporization of surface particles by microimpacts. 
It is noteworthy to recall that lunar ray craters (e.g. 
Tycho and Copernicus) show a moderate reddening even if 
their estimated age (or the corresponding exposure to the 
Sun) is relatively la rge, of the order of hundreds of Myr 
(jStoffler et al.l I2OOII ). Microimpact vaporization, however, 
might have nothing to do with the reddening of asteroids, 
perhaps due to much lower mean impact speeds in the main 
belt (5 km/s) compared to that on the Moon (18 km/s). 
Moreover, even if we decide to restrict our analysis to the 
asteroids, severe differences come out in the weathering 
properties of objects that should be -in principle- rather 
similar (e.g., Eros and Ida) . According to conclusions of 
iGaffevI 1I2OIOI ). it is not easy, and probably potentially 
misleading, to define a simple -uniparametric- asteroid 
weathering scenario. 

On the other hand, both the solar wind and micrometeorite 
bombardment weathering processes coexist for all asteroids. 



even if their relative efficiency may vary depending on 
the heliocentric distance. The different effects of space 
weathering should be foremost dependent on the different 
properties of the impacted surfaces (composition, texture, 
presence of regolith and so on). Unless we are able to 
find a well-defined selection criterion, which divides the 
asteroids into different groups, with different values of all 
the variables affecting the individual reaction of the bodies 
to the space environment, the weathering properties of 
individual bodies are difficult to assess, but the possibility of 
a statistical analysis and the search for a unitary statistical 
model remains open. 

In this series of papers we discuss the role of mini- and 
micro-collisions, as a potential solution to the two-timescales 
conundr um (as already sug gested, at a very qualitative 
level, in IPaolicchi et al.ll2009l ). 

The role of gardening in affecting weathering properties 
has been discussed in the literature (see, for instance, 
iGil-Hutton et al.ll2002l), and a quantitative model has been 
introduced by IWillman et al.l (|2008l . l201(il ). However, their 
models, based on some ad hoc assumptions (see later) 
seem to not be capable of fitting together the properties 
of the asteroids with those of ordinary chondrites ( PCs). 
In a recent improvement l|Willman fc Jedickd 1201 ll ). the 
fit with ordinary chondrites could be attained, but the 
reddened surfaces of young asteroid families do not fit 
with the inferred space weathering timescale. Note that, 
however, not all S-type asteroids are expected to be 
genetically linked to ordinary chondrites. Indeed, several 
observed S-types have inf erred compositions incompatible 
with ordinary chondrites (|Gaffev et al.l Il993l : iMarchi et al.l 
i2005i^ : while laboratory analysis show that OCs may have 
originat ed from a comparat ively small number of parent 
bodies (|Burbine et al.l |2002| ). Therefore, we caution that 
the average unweathered spectral slope of S-types may be 
different from that of ordinary chondrites. 
Apart from the technical details, the essence of the problem 
can be qualitatively sketched in the following way. We 
have, on one side, some evidence for a short (~ 1 Myr) 
timescale for substantial reddening: experiments on ion 
implantation and the observed reddening of some of the 
youngest asteroid family members. On the other side, we 
have also hints for a longer timescale (~ 10^ — 10"^ Myr), 
both from micrometeorite bombardment experiments and 
from a significant continuation of the reddening on old 
asteroids. However, at a microphysical scale, both processes 
must be present simultaneously. Thus if the short timescale 
holds, there is no confiict, while, in the opposite case, one 
should understand why the faster process does not work. 
Moreover, if we introduce a de-weathering effect due to 
collisions, we can have two extreme cases: if the collisional 
timescale is by far shorter than that of the weathering, 
no significant reddening can happen at all; if it is by far 
longer, the collisions do not much affect the reddening. 
Thus, the role of collisions, within this scheme, may be 
relevant only if the gardening has a timescale comparable 
to the weather i ng; th is may be the case, as suggested by 
I Willman etahl (|2010l ). but this is not enough to solve the 
two-timescales conundrum. 

Thus, the solution may require an additional effect to 
be considered. In our model, this effect is connected to 
the saturation of the spectral properties, which can take 
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place even in the presence of frequent gardening of the 
surface, due to the reaccumulation of the ejecta created 
by an impact; this reaccumulation can partially cover the 
surface with already reddened layers. The relevance of the 
reaccumulation may differ in importance, depending on the 
target's size: thus a different behavior for small vs. large 
objects has to be expected. 

In this paper, we show new observational results in support 
of the above mechanism. Detailed numerical simulations of 
the regolith evolution and impact processes will be deferred 
to the next paper in the series. 



2 CONFLICTING TIMESCALES 

2.1 The "short-time" camp and the problem of 
near-Earth asteroids 

The most explicit statement of a fast reddening of S- 
complex asteroids is presented in P2, where the analysis is 
restricted to a set of MBAs, namely members of dynamical 
families. The reddening timescale is obtained, in P2, by the 
correlation of the spectral slope of asteroids (corrected to 
eliminate compositional effects) and their family age. The 
timescale obtained, for most of the reddening, using data 
concerning two young MBA families (namely. Datura and 
Lucascavin), is of the order of 1 Myr. This value is consis- 
tent with the estimates from ion b ombardment laboratory 
experiments (|Strazzulla et al.ll2005l ). therefore the authors 
of P2 claim that the reddening of young asteroid families is 
due to the solar wind. However, they find also a slower fur- 
ther increase of the spectral slope, which, according to the 
authors, might be due to different physical processes. Note 
that the observed initial fast reddening is based on the as- 
sumption that the unweathered slopes of the S-type Datura 
and Lucascavin families are the same as that observed in or- 
dinary chondrites. Otherwise, the resulting reddening time 
dependence may be different, even showing no intrinsic ev- 
idence of a double timescale. However, the existence of a 
fast reddening pr ocess and the observed peculiar proper ties 
of Q-type NEAs l|Marchi et al.ll2006bl : [Binzel et al.ll2010l ') do 
not support this extreme possibility. 

The reddening, if mainly caused by the Sun, has to be re- 
lated to the exposu re E to the solar wind, as defined in 
iMarchi et all (|2006al ): 



E '. 



a2(l 



For the data sample used in P2, the slope-age and 
slope--E plots are rather similar, since the range of helio- 
centric distances is rather narrow (2.23 — 2.87 AU). The 
slope--E relation, obt ained from the same data used by 
IVernazza et al. I (|2009l l. is represented in Fig. [1] On the same 
plot, we a lso report the a v erage data point from the NEAs 
dataset of IPaolicchi et aP (|2007l l. For a more detailed com- 
parison between NEAs and MBAs, we correct the average 
NEA's data point for surface composition and perihelion 
de-reddening effects due to close encounters with terrestrial 
planets. For the composition effect, a mean relative abun- 
dance of olivine an d pyroxene ol/(ol + o px) ^ 0.7 was es- 
timated for NE As (IVernazza et all 120081 '). Thus, using the 
formula given in IVernazza et al.l (|2009l ) to correct the slope 



for a different composition, the mean spectral slope of NEAs 
has to be slightly increased by -|-0.07/im~^ (Fig- tt])- This 
data point was further corrected in order to take into ac- 
co unt the periheli o n de- reddening correction, as introduced 
in IPaolicchi et al.1 (|2007l ). 

We see that, in all cases, the NEA data points are too low 
compared to MBAs for a similar exposure. On the other 
hand, the average NEA slope is similar to those of the 
Datura and Lucascavin families, thus implying a similar ex- 
posure. Due to the different heliocentric distance, this expo- 
sure corresponds to a typical age of the order of 0.1 Myr for 
NEAfl 

This timescale is very short compared to all the relevant 
evolutionary timescales concerning NEAs: 

(i) their average ages (or collisional lifetimes) taking also 
into account the time passed in the main belt (the lat- 
ter may exceed the dynamical lifetime of a typical NEA 
and have a dominan t role in determining the reddeding; see 
iMarchi et "al]l2006al '); 

(ii) their typical lifetimes as NEAs; 

(iii) the time i nterval between significant de-weather ing 
close encounters (|Marchi et al.ll2006bl : fBinzel et al.ll2010l) . 

In other words, the NEAs should be almost completely 
reddened, a conclusion that is in stark contradiction with 
observations. 

A drastic solution to the problem assumes that the short 
timescale for reddening estimated by P2 is a fluke, and the 
correct timescale is of the order of hundreds of Myr. In 
this case, the NEAs problem would be solved; however, to 
match ordinary chondrites and young family objects is far 
less easy. In fact, the typical slope of the youngest known 
asteroids is, for the most part, significantly redder than 
that of the ordinary chondrites. 

It is also possible to suggest that the timescale estimated 
in P2, and based only on a couple of observational points, 
is underestimated by, say, one order of magnitude. In this 
case the de- weatherin g effect due to close planetary passes 
(|Marchi et al.l l2006bl : iBinzel et all 120101 ) might do the 
job. This possibility should be enforced by weakening, as 
discussed before, the link between the slopes of ordinary 
chondrites and unweathered S- types. 

Another, and more intriguing, possibility, will be discussed 
in greater detail below, namely that the main difference 
between NEAs and MBAs should be the size: most NEAs 
are small and, if the gardening refreshing is rather fast but 
can be significantly weakened by the self-reaccumulation of 
the ejecta, then small bodies cannot become too red. See 
Sect. 3 for further discussion. 



2.2 The "long-time" camp 

The long-tim e camp is mainly represented bv lWillman et ah! 
(200l, I2OIOI ). Their timescales are of the order of several 
hundreds of Myr, and a different laboratory counterpart 



^ A similar conclusion is obtained by also rescaling the differ- 
ence in exposure between the young MBA families and a typical 
NEA (a ~ 1.5 AU; e ~ 0.4). The reddening timescale of 1 Myr 
suggested by P2 becomes 0.1 — 0.2 Myr for a typical NEA. 
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(nam ely, micrometeorite impacts) is suggested l|Hiroi et alj 
The long-time camp is partially supported also by 
the results of our group (PI). We have found that a residual 
reddening (20%) takes place over very long ages, even of 
the order of 1 Gyr. However, in our model most of the 
reddening takes place in the fi rst period, according to the 
analysis and plots presented in iPaolicchi et~al] l|2007l ). 
Moreover, there is evidence for a fast reddening of some 
young objects (such as those of the Karin family), in 
substantia l agreement w i th th e short-time camp sug- 
gestions l|Paolicchi et alJ |2007| ). Finally, the apparent 
Sun- dependence of the weathering (the slope is more 
strongly dep e ndent on the exposure than on the age: 
iMarchi et all l|2006al )') strongly supports the dominance 
of the solar wind-driven ion-implantation processes, even 
if it cannot be considered as final, unequivocal proof: in 
fact, the micrometeoritic bo mbardment m ay also depend 
on the distance from the Sun (|Cintalalll992l ). as claimed by 
IWillman etall (|2010l ). Note, also, that the total fluxes used 
in experiments are tuned to those expected to come from 
the Sun, even if the laboratory rate of ion bombardment 
is, obviously, by far larger. Thus, if the ion-implantation 
process is not the main cause of weathering, one should find 
a theoretical explanation of it. 

In summary, we are facing a process characterized by two 
reddening timescales: the former, of the order of a few Myr 
(or even less), is characterized by significant reddening of 
several young (and typically small) objects; the latter, of 
the order of lO'^ — 10^ Myr, is characterized by a further 
reddening, towards saturation of the effect. While there 
is no reason to suggest that two different microphysical 
effects are at work (why is the faster one not able to redden 
after the first few Myr?), the possibility of a complex 
process, in which weathering and de-weathering effects are 
simultaneously active, seems to deserve serious scrutiny. 



3 GARDENING VS. WEATHERING: THE 
MODEL 

The timescale of collisional gardening is not easy to com- 
pute, since it depends on several (partially unknown) pa- 
rameters. This is also the basic reason for the forthcoming 
numerical simulations (see Sect. 5). Recent estimates, pre- 
sented in the l iterature, differ by several orders of magnitude 
jWillman et a l. 2010; Mclita ot al. 20Q9;) and it is not easy 
to solve the apparent discrepancies. In this Section, rather 
than present a new computation, we identify the main pa- 
rameters of the problem. Essentially, the basic questions are: 

• Ql. For a given asteroid size, what is the ratio between 
the collisional disruption timescale (ted) and the global 
resurfacing timescale -not automatically entailing a general 
spectral refreshing (see later in the text)- due to an individ- 
ual impact (tgr)? 

• Q2. Given the size distribution of the projectiles, what 
is the gardening timescale (tga) due to all projectiles? The 
result is obtained by integrating between a maximum pro- 
jectile size (Dmax) and a minimum size (-Dmin). When is the 
integral dominated by the value of Dmin? 

• Q3. What is the minimum impactor size to be consid- 
ered? 



• Q4. If a crater is formed on a reddened region, what is 
the fraction of the ejecta which, retailing onto the asteroid 
surface, gives rise to a -partially or completely- reddened 
surface? This possibility has not been taken into account in 
previous computations. The relevance of the effect depends, 
obviously, on the size of the asteroid, thus introducing dif- 
ferences between small and large objects. 

L et us consider Ql first . Acc ording to the analys is presented 
in lWiUman et al.] (I2OI0I ') and lMelita et al.l l|2009l '). the size of 
the resurfaced region (drcs) is roughly proportional to that 
of the impacting projectile (D), namely dros — ceD. Thus the 
resurfaced area is times larger than the projectile cross 
section. The value of a depends on various parameters, but 
the linear relation seems rather reasonable and robust for 
the typical conditions of asteroidal impacts: in fact the size 
of the resurfaced re gion is proportional to the size of the 
crater l|Meloshlll989l ). and this latter is proportional to that 
of th e projectile (at least in the strength regime; iMeloshI 
Thus we can adopt this as a basic rule. Consequently, 
it is possible also to discuss the first question: a global resur- 
facing follows from the impact of a projectile whose size is 
about 2/ a times that of the target (the factor 2 comes out 
from the factor 4 relating the surface of a sphere and the 
area of a circle with the same radius). 

The value of a is not easy to estimate. According to 
IWillman et~aLl (|2010l ). and references therein, the size of the 
crater dcra — 137?, and the size of the resurfaced region is 
about 2.3 times larger than the crater size. The mean size 
of the resurfaced region is thus about 30 times that of the 
projectile (or, in terms of area vs. the cross-section of the 
projectile, we have a factor around 10'^). However, the size 
of the crater might be larger: for instance, the scaling laws 
suggest a value approxima tely twice as large: mor eover, the 
Deep Impact experiment l|Richardson et al1l2007^ suggests 
a transient crater, say, more than 100 times larger than the 
projectile. On the other hand, to compute the ratio between 
the resurfaced area and the crater are a is even more difficult : 
the estimates might be a fa ctor ~ 5 (|Willman et al.|[2010l ). 
or ~ 25 (.Melita et al.ll2009l). which is intermediate between 
those suggested by Gil-Hutton et all l|2002l ). according to 
which the ratio between the sizes (to be squared to con- 
vert to areas) is between 2 and 10. The above uncertainties, 
all together, may affect the value of a, increasing it from a 
minimum value of about 30 (see above) , by even more than 
an order of magnitude, consequently strongly decreasing the 
gardening timescale. 

Therefore, the minimum projectile/target mass ratio caus- 
ing a complete gardening is ~ 3 x 10~^ in the conservative 
case, while going down to 10~* in the extreme opposite; 
the former value is , according to cu rrent collisional theories 
(iBottke et al ]|2005l : lHolsappldl2009D , smaller than that caus- 
ing a catastrophic disruption but lar ger than that causin g 
a complete shattering of the target (jWillman et al.|[2010l ). 
If, as usual, the size distribution is a monotonic decreasing 
function of size, it entails tgr < ted (maybe igr <C ted) - More- 
over, especially when considering large targets, the role of 
gravity has to be taken into account. 

Note that in order to evaluate the global effects of colli- 
sions, it may also be important to consider the so-called 
"global-jolt", namely regolith displacement over the entire 
target surface due to impact-generated seismic waves. These 
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effects can be estimated according to lO'Brien et al] l|2006h . 
by the following equation: 

i3*(<7/pO«KV«')''" (1) 

where D* is the size of the projectile producing, as a con- 
sequence of the impact, an acceleration g (in units of the 
acceleration of gravity g^) to all surface particles, and is a 
function of the target diameter dt and of the impact speed 
V. Assuming that similar consequences follow from a similar 
value of g/gt, the above formula indicates that the size of a 

5/3 

projectile able to produce global-jolt scales as . There- 
fore, for the same impact speed, smaller asteroids are af- 
fected by progressively smaller impactors than larger ones. 
Previous arguments showed that the physics of the impact 
processes plays a major role in understanding the collisional 
disruption rate and the rate of global resurfacing. Since the 
outcome of a collision depends on the physical properties of 
the target, maybe a single recipe for all asteroids does not 
hold. More likely, small rubble-pile asteroids will have a dif- 
ferent response compared to large rubble piles or monoliths. 
The answer to question Q2 requires the choice of an aster- 
oid size-frequency distribution. The problem is not simple 
and will not be discussed in detail here. However, we can 
take a simple power law, such as the traditional "Dohnanyi 
slope" (Dohnanyi, 1969) obtained by simplified assumptions 
regarding a stationary collisional cascade: 

dN = Am'" dm = A' D'"' dD (2) 

where A and A' axe constant quantities and m is the mass. 
In general q' = 3q — 2; for the Donhanyi slope q = 11/6, 
q' = 7/2. If we combine this assumption with the above- 
quoted ansatz dres = aD (and thus Stcb ~ {'k / 4)o? \ Srcs 
is the resurfaced surface) we obtain the relation: 

dS.es = S',es(-D)diV oc 0'^"+'^ dD. (3) 

The integral has to be performed between a maximum 
size Dmax of the order of that causing the collisional dis- 
ruption (with a corresponding timescale tcoii), or, better, 
of the order of the projectile size that causes global resur- 
facing Dgi, and a minimum Dmin to be defined. With the 
above assumptions, we obtain something proportional to 
^min^ ' — -D^/'"'"^. If so, whenever q > 5/3 (or, equivalently, 
whenever the differential size distribution has an exponent 
steeper than —3), the small impacts dominate the timescale: 

tgs. — tgr{Dgi-/Dn^in)'^" ^ < tcoll (-Dmax/-Dmin )^'' ^. (4) 

Assuming a Donhanyi slope (3g — 5 = 1/2), for an as- 
teroid of size a few km, which should be disrupted by a ~ 
km-sized projectile, the gardening time is of the order of 
tgr/20 (or less) if the minimum useful size of the impactor, 
to cause damage, is assumed of the order of a few meters (as 
in lwillman c t al. 2010), whi le decreasing below 10~^tcoii as- 
suming Dmin — 1 mm (as in lMelita et al.| [ 200!^ ). The values 
can change also assuming a different size distribution or al- 
tering other parameters of the model. 

Figure 2 compares the analytic distribution wit h that of 
MB A s derived by a collisional evolution model (Bottke et al.l 
[iooi). The plot also shows the MBA size distribution 
deriv ed by SUBARU observations (|Yoshida fc Nakamural 
l200'iD . which is valid down to ~ 1 km. All distributions are 
normalized at 5 km. The plot shows that the Dohnanyi slope 



is a good approximation of the MBA size-frequency distri- 
bution for the purpose of this work, although, there may be 
significant local slope variations. 

The discussion above shows that gardening is usually dom- 
inated by small impacts and that the estimate of the min- 
imum impactor size causing resurfacing (Q3) is crucial for 
estimating the gardening timescale. 

The last question of the list (Q4), however, may affect the 
above conclusions. In the case of an impact followed by the 
recapture of all the ejected fragments, one may imagine that 
a non-negligible fraction of the fragments (cr) will show an 
already-reddened surface. This fraction should be particu- 
larly large in the limit of a very small crater. We assume 
that a newly formed crater excavates material from pre- 
existing ejecta layers. This material will be mixed-up and 
then spread around to form a new ejecta blanket. Let /3 
be the volume ratio of reddened particles to non-reddened 
particles present in the ejected material (/3 = 0, 1 for non- 
reddened and totally reddened particles, respectively). Thus, 
the resulting fraction of reddened particles exposed in the 
new ejecta blanket is a = 0.5/3, assuming that half of the 
particles should fall showing the reddened surface. 
If we neglect this effect, the maximum attainable spectral 
slope, due to the combined effect of weathering (with a 
timescal e of tsw) and collisions, s hould be a function of 
tsw/tga (jWillman et al."2008'. 2010). A similar analytic toy 
model is discussed in Paolicchi ct al. (2009,). Conversely, if 
we take it into account, one can imagine a progressive slope 
saturation, with a timescale of the order of iga/o". However, 
in the meantime, larger impacts take place, affecting deeper 
-and fresher- regions; the resurfacing causes a partial refresh- 
ing (the fraction is related to the fraction of deep, unweath- 
ered ejecta) and the slope does not completely saturate. Af- 
ter a longer time, these layers also become completely or ex- 
tensively reddened, and the saturation is limited only from 
the consequences of even larger -and rarer- events. Thus the 
slope saturation may take an asymptotic behavior; one can 
guess that the final timescale to approach the slope s atu- 
ration may be of the order of ted (jPaolicchi et al.l 120091 ). In 
principle one might suggest that the reddest asteroids should 
be those whose typical age (or lifetime) is close to the age 
of the Solar System. Smaller asteroids should be bluer -for 
the reasons discussed above- than intermediate sized bodies; 
the same holds true for the larger bodies being farther from 
the asymptotic saturation. In the next Section we compare 
our ideas with observations. 



4 COMPARISON WITH THE OBSERVATIONS 

The discussion presented in the previous Section is rather 
challenging. Some suggestions may be verified only with 
the aid of a detailed model of surface evolution, requir- 
ing numerical simulations. In this Section we introduce 
some observational data that seem to support, at a semi- 
quantitative level, our ideas. 

One of the most interesting points in our model is, 
certainly, the different behavior between small and large 
asteroids. For the former, gardening is more effective, since 
the reaccumulation of partially reddened ejecta is absent 
or strongly reduced. If the asteroids are a little bit larger. 
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gardening is less effective, and the surface, for a given 
exposure, might be redder. This trend should continue as 
far as the asteroid size increases, reaching a plateau when 
the reaccumulation becomes massive or almost total. The 
value of the corresponding size depends on the physics of 
coUisional processes; it has certainly to be within the 1 — 10s 
of km size range, since reaccumulation requires the typ- 
ical speeds of the ejecta to be smaller than the escape speed. 

For larger bodies, as a first approximation, we imag- 
ine the spectral slope to remain close to the plateau value, 
but, if the arguments presented at the end of the last sec- 
tion are valid, the complete "slope saturation" is reached 
in a time which is of the order of the coUisional lifetime. 
Thus, if the lifetime is larger than the possible age of the 
body (which cannot exceed the age of the Solar System), 
the slope might be slightly under the saturation value. The 
transition is for bodies o f the order of a few tens of kilo- 
meters (see, for instance. iBottke et al.ll2005h . In Fig. [3] we 
show the exposure-corrected slope as a function of size for 
MBAs. The exposure-corrected slope is computed multiply- 
ing the observed slope of a particular object by the ratio 
between its estimated exposure (function of the orbital pa- 
rameters and of the age) and the mean (among all asteroids) 
exposure. Thus, the obtained exposure-corrected slopes are 
scattered around an average value with no correlation with 
the asteroid sizes. We find a statistically significant increase 
up to around 15 km (2-tailed probability of 2.2 x 10~^^), 
followed by a significant decrease beyond 15 km (2-tailed 
probability of 3.5 x 10~^). Note that the exposure, used to 
obtain the data represented in the figure, assumes constant 
orbital element s (a, e) and an age com puted according to the 
prescriptions of lMarchi et al.l (|2006al ) (applied to all objects 
regardless their family membership). However, for family as- 
teroids a different age estimate, obtained from the analysis of 
the overall properties of the related family, can be obtained. 
We verified that the quality of the trend is not affected by 
the definition of the age. With all the uncertainties in the 
model. Fig. [3] supports our ideas. 



5 NUMERICAL SIMULATIONS 

As discussed in the previous sections, the basic problem 
to be solved for a reliable assessment (or falsification) of 
the ideas sketched above, and partially supported by obser- 
vational evidence, is to understand what happens when a 
crater is formed on the surface of an asteroid. In particular, 
the critical issues are: 

• II. What happens when the projectile impacts onto a 
fractured (regolith, imbricated) surface and the size of the 
projectile is smaller or comparable to that of the surface 
components? (This projectile size is represented by Dmin in 
Sect. 3.) 

• 12. What is the size and velocity distribution of the 
ejecta (and, possibly, the size- velocity relation); in other 
words, given the size of the target, how many -and which- 
ejecta will fall again onto the target, and where (within the 
crater, close to the crater, everywhere on the surface)? 

• 13. Is there any effect due to the rotational properties 
of the ejecta? In other words, do they, re-falling, show ap- 
proximately the same external surface in about 50% of the 



cases (consistent with common sense), or is there some sub- 
tle reason for a different value? 

The handling of these problems, with the use of numer- 
ical simulations, might be successful, as indirectly shown 
by the possibility of explaining some observed feat ures of 
Eros in terms of a dynamical model of crater ejecta (|Durdal 
l2009h . 

Direct simulation of regolith dynamics is becoming more 
feasible with advances in computer sp eed and numeri- 
cal alg orithms. Sanchez et al. ( 2oTll) and iRichardson et ahl 
(|201ll ) have begun development of discrete element methods 
tailored for the low-gravity environment of small asteroids. 
Key advances include efficient collision handling, proper 
accounting of surface friction, and implementation of weak 
non-gravitational forces that may play a critical role in the 
evolution of surface regolith. 

Many parameters are uncertain, even if critical information 
and constraints may come out from the analysis o f space 
missions (see, for instance. IRichardson et ahl |2007| ). How- 
ever, new interest in sample return from and possible human 
exploration of small asteroids is spurring development of 
laboratory and computer experiments that may ultimately 
shed light on the weathering processes being discussed here. 
We will devote a forthcoming paper to the implementation 
of the required numerical simulations and to the discussion 
of the results. Simulations being developed will allow a 
portion of the granular surface of an asteroid to be modeled 
as a collection of discrete, possibly non-spherical particles 
in resting contact. Impacts and/or seismic shaking of the 
region will be simulated to determine the extent of ejecta 
redistribution and overturning as a function of impactor 
size, speed, and incidence angle. 



6 CONCLUSIONS AND FUTURE WORK 

The purpose of the present series of papers is to es- 
tablish a complete model of the space weathering pro- 
cesses for S-complex asteroids, in terms of a balance be- 
tween the reddening due to weathering and the refresh- 
ment of the surface due to coUisional processes, introduc- 
ing also the possibility of a reduced refreshment arising 
from the reaccumulation of coUisional ejecta. In this pa- 
per we have outlined the general features of the scenario 
and the most relevant uncertainties of the theory. In spite 
of these uncertainties, we have suggested a way to over- 
come the apparent conundrum of the relevant timescales 
coming out from the data. The suggested explanation is 
also -at least qualitatively- supported by a particular anal- 
ysis of the data. However, the overall complexity of the im- 
pact/cratering/ejection/reaccumulation processes requires a 
more detailed analysis in terms of numerical simulations, 
which we intend to perform. 

The realistic goal of this future study is a more quantita- 
tive estimate of the involved timescales, allowing a better 
general model and a more detailed statistical analysis of the 
data. However, we guess that the detailed interactions be- 
tween external disturbances and the surface properties of 
any individual asteroid may lead to a wide spread of results. 
Their interpretation, and the related possibility of deriving 
some relevant surface properties for a given object, requires 
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very sophisticated modeling. In this sense, the verification or 
falsification of the suggestions presented in the last Section 
may be of some preliminary use. 

Finally, we remark that a reliable model of space weathering 
of asteroids is not only interesting per se, but also represents 
a fundamental step to understanding the overall surface evo- 
lution of asteroids, including collisions, thus cratering and 
erosion, formation of regolith layers, and so on. 
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Figure 1. A slope-exposure relation corrected for chemical com- 
position, based on the data of Vernazza et al. (2009). Average 
NEA datapoints are also included (triangles), represented with 
its nominal age as estimated in Marchi et al. (20D6a, 2006b) and 
Paolicchi et al. (2007). The three points correspond to the nom- 
inal mean slope, to the same data corrected for the perihelion 
term (labeled as "pes"), and also with a further compositional 
correction due to the different mean value of the olivine-pyroxene 
parameter between MBAs and NEAs (labeled as "pcs-|-ccs" ) . See 
text for further details. 
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